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SU(2) lattice gauge theory is investigated where the traces of the Wilson lines at any lattice point 
and along each direction is constrained to zero. Hence, each of the lattice configurations possesses 
a vanishing density of heavy (anti-) quarks. The results are compared with those of pure SU(2) 
gauge theory which can be interpreted as the grand canonical realization of the heavy quark theory 
where only the ensemble average of the heavy quark density vanishes. 

The static quark anti-quark potential of the constrained theory is obtained from (spatially smeared) 
Wilson loops at zero temperature. We find that the potential coincides with that of pure SU(2) 
gauge theory (without constraints). Hence, the familiar "running" of the lattice spacing with j3 is 
recovered. 
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1. The role of the Polyakov line 

1.1 Heavy quark physics 

If the quark determinant is calculated for an arbitrary gluon field in the heavy quark limit, the 
correction to the pure gluonic action density is given by [||] : 

jl— in —/J —in . 

SheavyM <* <?~ tlP(x) + tlP J (x) , 

where m is the heavy quark mass, T is the temperature, and /I is the quark chemical potential. P(x) 
is the Polyakov line, i.e., the product of all time-like links forming a line parallel to Euclidean time. 

Taking the derivative of the partition function with respect to the chemical potential provides 
the expectation value of the baryon density. At vanishing chemical potential, we find 

<p) - (3 trP(x)} . 

In the case of SU(2), this density vanishes because the trace of any SU(2) matrix is real. For 
SU (N > 2), (p) vanishes if the center symmetry is realized: If the center transformation 

P(x) ->P'(x) = zP{x) , z = exp|/^m| , -N <m<N , 

is realized, real part and imaginary parts of (trP(x)} vanish. While a particular lattice configuration 
makes a non-trivial contribution to trP(x), the ensemble average causes its expectation value to 
vanish (see figure |2| for an illustration for the case SU(2)). This is no longer the case if T exceeds 
the critical temperature T c for deconfinement: In the latter case, the above center symmetry is 
spontaneously broken. The ground state suggested by the lattice ensembles seems to be populated 
with heavy quarks of Af-ality m s (where m s specifies the realized center sector). 

Lattice configurations of a standard simulation of pure SU (N) gauge theory with periodic 
boundary conditions can be interpreted as part of a grandcanonical ensemble where, on average, 
the density of heavy baryons vanishes. By contrast, if we demand that 

3 trP(x) = 

for each configuration, we sample a canonical ensemble for which the baryon density is identical 
zero for each ensemble member. 

1.2 Topological sectors of perturbative vacua 

It was pointed out by M. Schaden in [jl]] that even the Yang-Mills groundstate can fall into 
topological different classes on a torus. Let us briefly review his arguments. 

A perturbative ground state of SU(N) lattice Yang-Mills theory with Wilson action is given by 
a link configuration with 

-trP MV (x) = 1 , Vx, v>n , (1.1) 
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Perturbative vacuum: 



Trace Polyakov line: 



U 



U 




u 



u 



U = exp( i Q n^t 3 ) 
U 

N, U 
U 
U 



U constant 

U i (topol.) different from U2 for: cos( N t 1) 4 cos( N t 62) 
Figure 1: Topological different perturbative sectors on a torus for a SU(2) gauge theory. 

where P liV {x) is the plaquette built from the links U^(x): 

P MV (x) = U^(x)U v (x + n)U J( n{x + v)Ul(x) . 

Obviously, the condition (|0|) is satisfied for all constant links 

U^x) = U =: exp{/0V}, Vx 5i u 

and their gauge equivalents (see figure [IJ). Let us consider the Wilson line P in time direction 
(Polyakov line) for a lattice with N t links in time direction. Since the trace of the Polyakov line 
is gauge invariant for a lattice gauge theory with periodic boundary conditions, the perturbative 
vacua, specified by U, can be classified by 



K(U) = ^trP = ^tiexp{iN t G a t"} . 
In particular, two candidates Ui and U2 are topological inequivalent if fc(t/i) 7^ k(U2)- 



(1.2) 



2. Constrained Yang-Mills theory 
2.1 Theory 

In the following, we will show results from a simulation of SU (2) and SU (3) lattice gauge 
theory on a 4-torus using lattice configurations with periodic boundary conditions. The lattice 
configurations are subjected to the constraints: 



trP(x) =0, VI 
3 trP(x) =0, VI 



for 51/ (2) 
for St/ (3) 



(2.1) 
(2.2) 



Results will be shown for temperatures below and above the deconfinement temperature. 



The results are interesting for the following reasons: (i) The results will show whether YM- 
theory with heavy quark fluctuations suppressed (i.e. constrained YM-theory) yields the same re- 
sults as pure Yang-Mills theory with periodic boundary conditions. Due to spontaneous center 
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symmetry breaking at high temperatures, differences are expected at least for the deconfinement 
phase, (ii) If constrained and standard YM-theory yield the same physics at least in the confinement 
phase, the range of applicability of Schaden's center symmetric perturbation theory [jl]] might be 
even extended to SU(2) gauge theory, (iii) If constrained YM-theory is confining, the confinement 
mechanism might well be in the reach of a semi-classical treatment of the partition function based 
upon background fields with particular holonomies such as calorons [J| |J]. 

To be specific, the partition function of the simulation is e.g. in the SU(2) case given by 

Z = Js>U fj5(trF(x)) exp{s Wl i}. (2.3) 

One might ask whether the constraints in the 8-f unctions affect the bulk properties of the lattice 
ensemble at all: N t link elements contribute to a given Polyakov line. Hence, 3N t real numbers 
(for SU(2)) must satisfy one constraint. This seems to suggest that the constraint is irrelevant in 
the infinite volume limit. Note, however, that N t — 1 links can be gauged to unity. In addition, the 
remaining link can be chosen diagonal (up to a gauge transformation). This implies that in lattice 
gauge theory, there is only one physical degree of freedom contained in the N t links of the Polyakov 



line. From this arguments, we expect a significant influence from the constraints in (2.3). 



2.2 Numerical setup 

Firstly, a local update algorithm on the basis of M. Creutz heat bath algorithm was imple- 
mented. There, a particular link was chosen for the update, and the update of this link was per- 
formed according to its probabilistic weight, properly taking into account the constraint. 

satisfying changing the trace of the Polyakov line. Seemingly, this overemphasizes the con- 
straints, and one might fear that the algorithm is not ergodic. 

For this reason, we secondly constructed an algorithm where all links, "participating" in one 
particular Polyakov line, are updated simultaneously. We found that both algorithms yield the 
same results (at least for the observables shown below). Hence, ergodicity is not a problem for the 
present algorithms. 



3. Numerical results 



In a first step, we have calculated the static quark anti-quark potential at zero temperatures by 
means of (spatially) smeared Wilson loops. Our findings for the constrained case are summarized 
in figure ^| The striking result is that, despite the constraints, the potential is still linear rising 
at large distances. In addition, the potential agrees for each j8 separately with the potential from 
standard simulations. This implies that the scaling, i.e., the dependence of the lattice spacing with 
j6, is unchanged by the constraints. In addition, the potential was obtained for different orientations 
of the quark anti-quark pair relative to the cubic lattice ((110), (110), (110) directions). Our result 
shows that breaking of rotational symmetry by the underlying lattice is small. 

Subsequently the question arises whether constrained YM-Mills theory possesses a decon- 
finement phase transition at high temperatures despite the constraints which set the Polyakov line 
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Figure 2: SU(2): Distribution of the Polyakov line expectation values for temperatures below and above the 
deconfinement one (left panel). Static quark anti-quark potential as function of the quark anti-quark distance 
for constrained YM-theory at zero temperature (right panel). 
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Figure 3: SU(2): the sum of energy density £ and pressure p as function of temperature. 



expectation value to zero. In order to answer this question, we have monitored the "black body 
radiation" of gluons as function of temperature. In particular, we have calculated 

£ -±P = mN ?(P s -P t ) 

where e is the energy density and p is the pressure, P s , P t are space-like and time-like plaquettes, 
respectively, and /(j8) is a scaling function [Q]. Simulations were performed at fixed /3 while tem- 
perature was varied by changing N t . Our result for standard SU(2) YM-theory is shown in figure || 
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Figure 4: Same as figure 3 for the case of SU(3): standard SU(3) (left), constrained SU(3) (right). 

and agrees with the well-known results of the Bielefeld group [^] : The onset of deconfinement is 
signaled by a sharp rise of e + p. Figure || also shows our result for constrained YM-Theory. We 
find also sharp rise of e + p at roughly the correct deconfinement temperature T c . In addition, e + p 
strongly peaks at T c and seem to approach the Stefan Boltzmann limit from above. Unfortunately, 
we observe a strong cutoff dependence of e + p above T c , which renders a physical interpretation 
of the results cumbersome. 



Let us finally study the case of a SU(3) gauge theory with constraint (22). In figure |], the result 
for £ + p of standard YM-theory (left panel) is compared with that obtained from simulations of 
the constrained theory. In the constrained case we find that e + p is monotonically decreasing with 
temperature. No signal of a deconfinement transition is observed. Moreover, it seems that black 
body radiation from deconfined gluons is present at all temperatures. For an interpretation of these 
findings, it is important to notice that SU(3) center symmetry is broken by the constraints. Also 
center symmetry is merely associated with quark confinement, there seems also no confinement of 
gluons if center symmetry is broken. 
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